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Within the Replica Symmetry Breaking (RSB) framework developed by M.Mezard and G.ParisiQ 
we investigate the occurrence of structural glass transitions in a model of fluid characterized by hard 
sphere repulsion together with short range attraction. This model is appropriate for the description 
of a class of colloidal suspensions. The transition line in the density-temperature plane displays 
a reentrant behavior, in agreement with Mode Coupling Theory (MCT), a dynamical approach 
based on the Mori-Zwanzig formalism. Quantitative differences are however found, together with 
the absence of the predicted glass-glass transition at high density. We also perform a systematic 
study of the pure hard sphere fluid in order to ascertain the accuracy of the adopted method and 
the convergence of the numerical procedure. 

PACS numbers: 64.70.Pf 



I. INTRODUCTION 

Understanding the physical conditions leading to dy- 
namical arrest in colloidal suspensions and protein solu- 
tions is still attracting considerable interest [2} . Colloidal 
systems are often modeled as simple fluids of (large) 
spherical particles interacting via the hard core repul- 
sion plus, possibly, a short range attraction. The latter 
contribution may have diverse physical origins, accord- 
ing to the specific system considered, like the depletion 
mechanism or steric effects. 

This simple model has been investigated by simulation 
methods 0] and liquid state theories [J . The key dimcn- 
sionless parameter governing the physics of the system 
is the ratio between the range of attraction (A) and the 
hard sphere diameter (d) which, in the following will be 
taken as the unit of lengths. The equilibrium phase di- 
agram is now well known in a wide range of A and the 
stability of the fluid phase as a function of A has been 
extensively discussed in the literature. If the range A 
is such that at high densities many particles lie within 
the respective ranges of attraction, the overall effect of 
the attractions is expected [j| to be just a shift in the 
ground state energy of the system, while the dynamical 
and thermodynamical properties are mostly determined 
by the hard cores of the particles. Conversely, if the range 
of the attractive well is short enough, small displacements 
of the particles can break one or more "energetic bonds" 
and cause significant variations in the energy of the mi- 
croscopic state. While the hard sphere fluid is a purely 
entropic system, in this case we have both energetic and 
entropic contribution to the free energy whose relative 
importance can be tuned by varying the temperature 
and/or the density. 

The possible occurrence of dynamical arrest in this 
model has been studied via Mode Coupling Theory 
(MCT) which indeed predicts the existence of two dis- 
tinct glassy states, called "attractive" and "repulsive" 
glass, separated in the phase diagram by a first-order 



transition line terminating in a critical point @, Q • This 
finding might have relevant consequences on the experi- 
mentally accessible systems and it is therefore important 
to investigate the problem by alternative techniques. 

The occurrence of a glassy states in statistical physics 
was studied mostly by use of the replica approach in 
the framework of spin glasses When considering a 
system with quenched disorder, it is natural to intro- 
duce the whole ensemble of possible realizations of the 
quenched variables. For a variety of models it turns out 
that the values of many quantities of physical interest 
equal, in the thermodynamic limit, their average over the 
quenched variables. Subsequent studies revealed that in 
order to produce glassy behavior, an externally imposed 
quenched disorder is not essential because the frustration 
necessary for slow dynamics and dynamical arrest can be 
self-generated by inter-particle interaction. 

Throughout the formalism, different replicas are in 
principle indistinguishable and the equations for a repli- 
cated system are symmetric with respect to permutations 
of replicas; yet, replica symmetry breaking (RSB) solu- 
tions have proved to be the appropriate solutions for the 
low-temperature phases of several models. For such solu- 
tions, the correlation functions between different replicas 
do depend on the replica indices. 

A general framework often used to account for the 
slow dynamics of glasses is based on the analysis of a 
free energy F, considered as a functional of the aver- 
aged density profile [3. The key feature is that at low 
temperatures F is expected to develop a multi-minima 
structure and, if the free energy barriers are high enough 
with respect to UbT, the configuration space becomes 
essentially disconnected, leading to a non-ergodic behav- 
ior. The basic assumption underlying the mean-field ap- 
proach presented in this paper is the modeling of the 
glassy non-equilibrium dynamics through such a free en- 
ergy landscape. Hence, since the whole formalism of equi- 
librium thermodynamics promptly applies, we can refer 
to a glassy state as a proper phase of the system. The dy- 
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namic crossover observed experimentally is correspond- 
ingly replaced by a sharply defined phase transition. 

In this context, the introduction of replicas is a pro- 
cedure intended to unveil a possible multi-minima struc- 
ture of the free energy: it is natural to expect that for two 
replicas in the same local minimum of F the correlation 
is larger than for replicas in different minima. Hence, the 
appearance of a non-trivial pattern of correlation among 
replicas (i.e. the breaking of the replica symmetry) is the 
signature of the complex free energy landscape character- 
istic of a glassy phase. In the following, by studying the 
correlation functions between different replicas of a fluid 
system, we confirm the existence of a RSB phase tran- 
sition. The comparison of the resulting phase diagram 
with experimental and numerical data on dynamical ar- 
rest in systems with short range attractive interactions 
strongly supports the interpretation of such a transition 
as an ideal glass transition. 

We point out that by "ideal" glass transition it is 
usually meant a transition defined by a thermodynamic 
singularity [10|, With the replica method, two critical 
densities are obtained: the lower density signals the ap- 
pearance of a complex free energy landscape with no 
thermodynamic singularity. At the higher density the 
complexity vanishes non-smoothly and this gives a true 
thermodynamic transition. In it was shown that the 
replica approach we are going to use is not suitable for 
the description of the glassy phase and leads to incorrect 
results for the thermodynamic transition. For the time 
being we do not exclude that a more accurate numerical 
study of the thermodynamic transition might yield dif- 
ferent results (see table in section ITTT1 and consider that 
the grid used in [lj had only N — 128 points), however in 
this paper we limit our attention to the first transition, 
commonly known as "dynamical" glass transition. 

We first review the basic equations derived in Ref. [lj 
then we perform a careful study of the hard sphere sys- 
tem, in order to test the accuracy of the method and 
to estimate the numerical uncertainty of our results. Fi- 
nally, we investigate the effects of a short range interac- 
tion and compare our results with MCT predictions. We 
anticipate that, although the transition line we find is 
rather similar to those of MCT clearly showing the pre- 
dicted reentrant behavior, we have no evidence in favor 
of glass-glass transition, within the adopted approxima- 
tions. 



II. THE REPLICA METHOD FOR 
STRUCTURAL GLASSES 

According to the previous discussion, in order to imple- 
ment the replica symmetry breaking scheme for a model 
of fluid, we need a formal expression for the free energy 
functional of a mixture of identical copies of the original 
system. In the following, the indices a and b identify dif- 
ferent replicas. Such a generalized free energy functional 
F depends on the two body interactions among parti- 



cles U ab (r) (both intra and inter-replicas) and on the two 
point correlation functions g a b(i") = h ab (r) + 1, whose 
entries are labeled by replica indices: ordinary (intra- 
replica) 2-point functions on the diagonal and inter- 
replica correlation functions as off-diagonal elements. 
The original physical system can be recovered by set- 
ting n — 1. The key property of the functional F is 
to attain its global minimum when the pair correlations 
9ab(r) assume the physical value corresponding to the 
given interaction U a b(r). RSB does occur when a non 
trivial minimum (i.e. a solution with finite correlations 
among different replicas) is present in the physical limit 
of uncoupled replicas U ab (r) — U{r)8 ab . For a homoge- 
neous system, the general structure of the free energy 
density functional is 

F[U,g] = Y / B y l J d 3 xg ab (x)U ba (x) ~TS[g] (1) 

a.b 

The internal energy contribution is exactly given by the 
first term of Q while the entropy density S is known 
to be a functional of the pair correlation functions alone, 
which can be written [l2| as a sum of p-particle contri- 
butions: S = J2 P s p- By use of the Gibbs-Bogoliubov 
inequality it is easy to prove that the functional F 
defined in satisfies the requirements previously men- 
tioned. Unfortunately, a closed expression for the en- 
tropic contribution S is not available and we have to re- 
sort to some approximation. Following Ref. Q we keep 
the exact two body term of the excess entropy S2 

n „ 

s 2 = -k B ^ -y^ / d 3 x [g a b(x) \ng ab (x) - h ab (x)} 

a,b 

while approximating the residual contribution AS = S — 
S2 in terms of the pair distribution function: 

AS =~k B ±J ^_TrL(ph(g)) (3) 

where the function L(y) is defined as 

V 2 

L(y) = -]n(l + y)+y-^ 

and h(q) represents the n x n matrix whose elements are 
the Fourier transforms of the correlation functions h ab (r). 
This form follows form the resummation of the infinite 
class of non crossing diagrams in the formal expansion of 
the entropy as a functional of h(r). Collecting all terms 
together we find the final expression for the free energy 
functional which, in the physical limit of uncoupled repli- 
cas becomes: 

« n 

2/3F = p 2 j d 3 x^ (g ab {x)[\ng ab {x) + f3U(x)S ab ] 

a , b 

~ Kb{x) ) + J ^ TlL{ph{q)) (4) 
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FIG. 1: Replica symmetric (left) and 1-step replica symmetry 
breaking (center) form of the correlation matrix for a repli- 
cated system. The black squares represent single elements 
(<?*), while the dimension of the n x n matrix (light-gray for 
go) and of the m x m blocks (dark-gray for gi) is, in general, 
not fixed. In the figure, n — 6 and m = 3. Since in the physi- 
cal case of zero replica coupling the correlators corresponding 
to the light gray area vanish, the numerical computations can 
be simplified by considering only one of the blocks (right). 



In the case of a single component {n = 1), the min- 
imization of the free energy functional J1J leads to the 
usual HNC integral equation for the pair correlation: 



]ng(x) 
W(x) 



PU(x) 
d 3 p 



W(x) 



; Ph 2 (p) 

1 + phip) 



(•5) 



Yet, by keeping n unspecified, a second equation is ob- 
tained: the appearance of multiple solutions for this new 
equation signals a phase transition. Given the replicated 
free energy, the fundamental ansatz is a 1-step RSB. As 
we said, in general, a RSB directly mirrors a change in 
the "effective topology" of the configuration space which 
takes on a multi-valley structure, with "disconnected" 



single-valley domains. By choosing a 1-step RSB we im- 
pose only two possible patterns for the functional form 
of the correlations between replicas: either two replicas 
are in the same minimum of the free energy or they are 
in different minima; in the latter case we assume that 
the correlation of every pair of replicas is the same, no 
matter what the two minima are. In particular, for a 
structurally disordered material, we expect that any spa- 
tial correlation between two replicas in different minima 
should vanish: this conjecture follows from the obser- 
vation that the correlation between different disordered 
configurations, despite possible local similarities, should 
vanish in the thermodynamic limit, when averaged over 
the infinite volume of the system. A formal argument for 
the relevance of the 1-step RSB in structural glasses and 
a physical interpretation of m (defined below) in terms 
of effective temperature is presented in |l4j . In terms of 
correlation matrices, the standard procedure is to group 
the n replicas into n/m sets with m replicas each. Ideally, 
the m replicas within a block are in the same minimum 
of F. The "topological" structure of the matrix of corre- 
lation functions g a b is pictorially depicted in Fig. 1: The 
black squares represent the ordinary two-point correla- 
tion functions for a fluid (5*); the correlation functions in 
the light-gray blocks (50) are assumed to be unity (differ- 
ent minima are unrelated) and the correlation functions 
in the dark-gray blocks (gi) may be non trivial (repli- 
cas in the same minimum). Since each term in Q has 
a block-diagonal form and ho = 0, the free energy per 
replica turns out to be independent of n and the prob- 
lem reduces to the study of any of the n/m blocks of 
dimension m x m |l5j (figEt). Expressing F in terms of 
and g\\ 



2/3F 
np 2 



d 3 x( <7*(ai)[ln g*(x) + (3U(x)] — h*(x) — (1 — m)[g 1 (x) \ngi(x) — h\(x)] 



+ ph*(q) - (1 - m)ph 1 (q)} 



h*{q) hl(q) _ 1 - to 
p 2 mp 2 



d 3 ( 

(27r) 3 \ mp 2 

hi(q) 



ln[l 



ln[l + pK(q)-phl(q)}-(l-m)^f>) 



(6) 



The stationary condition with respect to is the usual 
HNC equation J5j, while, due to a simplification of fac- 
tors (1 — to), the extremal equation for the inter-replica 
correlation function gi is non-trivial also in the "physi- 
cal" case to — > 1: 

In Si (a;) = W 1 (x) 
Wl ( q ) = p[h*{q) - M<?)] 2 ^ 

1 + pK (q) 1 + p[h* (q) - hi (q)] 

Because of the non linear structure of equation l(7|l. we 
can now hope to obtain non-trivial solutions g\ (x) which 
would otherwise be ruled out by the replica-symmetry 
of the partition function. Indeed, we will show that the 



extremal equations for the free energy admit non- 
trivial correlations between replicas (i.e. a non trivial 
hi(x)), even for vanishing inter-replica coupling. The 
typical behavior of such a solution is shown in figure 2. 

To summarize, starting with a one component model 
of fluid, by introducing m identical copies and taking the 
limit m — > 1, after simple algebraic manipulations of the 
approximate free energy functional 10}, we obtain the 
usual HNC equation for the two-point function and one 
additional equation which describes correlations between 
different replicas. Since for to = 1 an inter-replica corre- 
lation has apparently no meaning, we could naively inter- 
pret h\(x) as an infinite time limit of a self-correlation 
function: The presence of a non trivial solution h\{x) 
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suit, we proceed as follows: starting from the "liquid" 
solution at low density {p < p t ), we introduce a finite at- 
tractive coupling e among replicas 0] and subsequently 
increase gradually the density until the continuous bifur- 
cation takes place. We then follow the glassy solution 
and finally reduce to the coupling between replicas, ob- 
taining a non-trivial correlation function h\(x) even for 
vanishing inter-replica coupling. The transition density 
Pt can then be pinpointed by decreasing (at zero inter- 
replica coupling) the density until the glassy minimum of 
the free energy disappears. The values we report in the 
following for the transition densities are the first ones 
for which our algorithm does not find a solution. The- 
oretically, the introduction of a coupling among replicas 
provides a way to discover multiple degenerate states, in 
this case it is motivated only by the computational con- 
venience to obtain a continuous bifurcation of the new 
minimum of the free energy from the known one. 



FIG. 2: An example of correlation functions for sticky spheres 
close to the dynamical glass transition. The thermodynamic 
state is p = 1.1878 T — 0.6 and the data refer to the system 
with A = 1/64. 

points to ergodicity breaking. This point is further dis- 
cussed in section flVAl 

A. Solution of the equations 

For any density, hi(x) = is a solution of equation 10. 
By linearizing (7J around hi(x) = it can be shown that, 
if another solution exists, it does not bifurcate from the 
trivial one. However, if we introduce a strong enough [l6| 
coupling between replicas, the linearized equation does 
not forbid a continuous bifurcation of solutions and a sec- 
ond solution has indeed been found. When the coupling 
is eventually turned off, one of these solutions converges 
to h\{x) = ("liquid" solution), while the other one 
maintains a spatial structure ("glassy" solution). This 
signals a multi-minima free energy, which we interpret 
as a non-ergodic phase. In fact, the free energy is now 
a functional of and g±: the occurrence of a second 
minimum with g\ ^ 0, mirrors the appearance of many 
minima in F[{p(x))], |17|. The glassy solution of equa- 
tion Q) does not exist for all values of the control pa- 
rameters but only in a region of the density-temperature 
plane, and the boundary between the 1-solution and the 
2-solutions domain in the phase diagram marks the glass 
transition. 

The two-point functions h^(x) and h\{x) are computed 
by the numerical solution of equation Q) . We first intro- 
duce a mesh in the radial coordinate, thereby reducing 
the unknown function h\ (x) to a set of N discrete values 
{h\(x) I— > hi(xi) i = 1, .., N). Then we use the iterative 
Newton-Raphson method for solving the resulting set of 
non linear equations . In order to find a non-trivial re- 



III. THE HARD SPHERE FLUID 

The glass transition in the hard sphere fluid has al- 
read y b een studied by the replica method in Refs. 
and [l6( . Here we consider the dependence of the criti- 
cal density on the choice of the discretization parameters 
used in the a numerical solution, and we provide an esti- 
mate of the transition density p t . 

Given the number of points N for the discretization of 
the correlation functions, the only arbitrary choice is the 
mesh size a = 1/Nd (where Nd is the number of points 
used to describe a hard sphere diameter) . It then follows 
that we can describe a correlation function h(r) only for 
r G [Q,N/Nd]. The physical limit should require both 
a -> and N/N d -> oo. 

After choosing a value for N d , we proceed by calculat- 
ing the transition density for grids with more and more 
points N while keeping N d fixed: for N large enough, 
the value of the transition density becomes independent 
of N and we can label it as pt(Nd). Repeating the pro- 
cedure for larger and larger A^'s the sequence of values 
we obtain converges to a finite value. Our results for the 
hard sphere fluid are summarized in table [I] 

It is clear that only when N/Nd ^ 4 the values pt(Nd) 
are almost independent of N. If we now consider the 
values obtained for the 8192-point grid, by arbitrarily 
interpolating the values pt{Nd) with a function of the 
form: 

Pt (N d = 2 a )=p t -Aexp(-ka) (8) 

we can extrapolate a limiting value pt- 

Depending on the inclusion of the first 5 or 6 points in 
the interpolation (a = 5, .., 9 or a = 5, .., 10 respectively, 
the boldface entries in the table), we obtain p t = 1.1699 
or p t = 1.1695 (corresponding to a packing fraction 
r\t ~ 0.612); this should be considered the final predic- 
tion for the transition density of a hard sphere fluid in 
HNC approximation. This result is rather close to the 
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1.1560 


8 


256 






1.1553 


1.1753 


1.1627 


1.1636 
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1.1578 
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1.1804 


1.1680 
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2048 












1.1595 


1.1812 
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4096 














1.1598 



TABLE I: Glass transition density pt for different choices of 
the discretization parameters: total number of mesh points 
N and number of points inside the hard sphere diameter Nd- 



1.18 



1.16 




1.14 



1.12 



Log 8 N d 

FIG. 3: Dependence of the critical density on the grid: the 
dots represent the data in the last column of the table, the 
dashed line is the best fit using equation JHJ. 



random close packing limit {q rcp ~ 0.64) and is consid- 
erably higher than the prediction of MCT (?y t ~ 0.516) 
0. We note that recent experiments suggest that, in the 
absence of gravity, the glass transition in "hard sphere 
colloids" is in fact remarkably close to r\ rcv [l9j . 

For the sake of completeness we obtained a similar ta- 
ble for the thermodynamic glass transition also (up to 
a 4096-point grid): the extrapolation of the asymptotic 
density yields p g ~ 1.189, which is equivalent to a pack- 
ing fraction rj g ~ 0.623. 



IV. THE SQUARE WELL FLUID 

In order to model a colloidal suspension, we now inves- 
tigate the effects of the presence of a short range attrac- 
tion by adding an attractive square well to the repulsive 
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FIG. 4: Phase diagram for the square well fluid for different 
well widths A. The dashed line represents the hard sphere 
transition density computed with the same discretization pa- 
rameters. 



hard core. The interaction potential is thus: 



oo r < d 
U{r) = { -U d<r<d- 
d + A < r 



(9) 



As usual, we take d as the unit of lengths. The appli- 
cation of the replica method is straightforward, but the 
phenomenology is now potentially richer, since both the 
density and the temperature are meaningful thermody- 
namic parameters. By inserting the potential @ into Eq. 
(J5J and applying the procedure described in section Til Al 
for different temperatures, we determine a whole transi- 
tion line Pt(T) in the p-T plane. 

In figure 4 we summarize the numerical study of the 
short range limit of the square well fluid. The results 
agree with the qualitative picture following from the 
above considerations: at high temperatures the glass 
transition line in the p-T plane is almost vertical and, as 
for the hard sphere fluid, and the glass transition is gov- 
erned by excluded volume effects. At low temperatures, 
the transition line displays a strong dependence on the 
temperature and only a slight dependence on the density: 
the dynamical arrest is driven mainly by the stickiness of 
the spheres. The replica method correctly captures both 
regimes, extending the transition line to densities well 
below p t . Moreover, in a region of the phase diagram, 
the "competition" between energy and entropy stabilizes 
the liquid phase at densities larger than p t . This phe- 
nomenon, commonly referred to as "reentrant behavior" , 
is also reproduced. 

The computations were performed on a 4096-point grid 
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(512 points for a hard sphere diameter). Because of the 
very short range of the attraction, few points were avail- 
able to describe the narrowest wells. Due to the non lin- 
earity of the integral equations, an estimate of the errors 
introduced by our discretization is quite difficult; yet, the 
qualitative behavior of the transition line seems to be well 
established. As an empirical test of the reliability of our 
results for wells described with very few (2 or 4) points 
on the grid, we repeated the computations for part of the 
phase diagram in figure 4 using a looser grid (2048-point 
grid, 256 points for a hard sphere diameter and a cor- 
responding rescaling by a factor of 2 for the attractive 
wells). The numerical results agreed within 0.2%. 

The convergence of our algorithm becomes rather del- 
icate as we follow the liquid-glass transition line to low- 
temperatures/low-densities and we defer to future stud- 
ies the behavior of the transition line in such a limit. An 
intrinsic limitation of our approach, comes from the equi- 
librium character of the replica method, which cannot be 
applied below the spinodal decomposition line. 

As already found in the pure hard sphere fluid, the 
transition densities are significantly higher (~ 15%) than 
those obtained by MCT [y| and the corresponding critical 
temperatures are lower, roughly by a factor 2. The reen- 
trant behavior of the glass transition, found by MCT, 
is confirmed by RSB method when A is sufficiently 
small. However, while in MCT the reentrance appears 
for A < 5%, in our calculations a similar shape is fond 
only for A < 1.5%. 



A. The glass-glass transition 

An interesting phenomenon that might be within the 
scope of the replica approach is the glass-glass transition 
predicted by MCT. In particular, in a dynamical theory, 
the localization length of a particle can be defined: a — 
lim t _ >oc (|r(i) — r(0)| 2 }. In a solid or in a glass, a attains 
a finite value, while it diverges in the fluid phases. It has 
been shown [|| that a undergoes a discontinuous change 
across the glass-glass transition line, decreasing from its 
value in the repulsive glass to a smaller value comparable 
to the attractive well width, in the attraction induced 
glass. 

In our formalism, this information is not directly ac- 
cessible. However, a hint on the local structure of the 
glassy state can be obtained by defining the quantity 
P(r) = 4irr 2 gi(r) which represents the probability (per 
unit length) that two particles belonging to different 
replicas lie at distance r one another and can be thought 
of as the infinite-time limit of a dynamical self-correlation 
function |2fT |. In figure 5 we plot such a probability for 
values of the control parameters which, according to the 
MCT, should identify distinct glassy phases. In all cases, 
the probability distribution is uniformly spread over a 
region of the order of a hard sphere volume. Even con- 
sidering the local maximum of P(r) for r ~ 0.05 (see 
figure 5), comparable to MCT localization length a, its 
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FIG. 5: P(r) for the same density (p = 1.19) and A = 1/128 
but different temperatures. 
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FIG. 6: Detailed structure of P(r) at short distance for the 
data of Fig. 5. As the temperature decreases at constant 
density p — 1.19 (see the vertical cut in figure 4) no singular 
behavior of P(r) is observed. 



position varies smoothly with the control parameters and 
its value is always close to the one obtained for the pure 
hard sphere system. We cannot exclude that a different 
choice for the replica symmetry breaking scheme (e.g. a 
2-step RSB) might lead to a well-defined glass-glass tran- 
sition: the 2-step RSB would then represent the splitting 
of the free energy minima revealed with the 1-step RSB 
into sub-minima representing a further trapping of the 
caged particles within the range of the attractive wells. 
This possibility is currently under investigation. 
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V. CONCLUSIONS 

In this paper we first investigated the occurrence of 
RSB in a model of hard spheres, providing an accurate fi- 
nite size scaling for the transition density. The mean-field 
replica approach implies a sharp (dynamical) glass transi- 
tion; the asymptotic result we obtained for hard spheres, 
at a packing fraction of r]R ep i ~ 0.612 should be compared 
with the MCT value tjmct ~ 0.52 and the commonly 
accepted value rj ~ 0.58 for the slowing down of the dy- 
namics of colloidal suspension and simulated hard sphere 
systems (even though experiments conducted in micro- 
gravity conditions 19] have recently raised this value, 
significantly closer to the random close packing limit). 
By introducing attractive interactions we obtained the 
glass transition line in the density-temperature plane. 
The expected reentrant behavior related to the change 
between "repulsive" and "attractive" glass is reproduced 
within the RSB approach. Quantitative discrepancies 



with respect to MCF are found for the square well fluid, 
analogously to the hard sphere case. Moreover, contrary 
to MCT, the attractive and repulsive glassy regimes are 
smoothly connected with no sign of a sharp transition. 
These results have been obtained within a simple ap- 
proximation of the more general RSB method: we just 
considered the one step RSB scheme and we evaluated 
the entropy functional by a HNC-like expression. It is 
clearly desirable to go beyond these limitations in order 
to check the qualitative and quantitative stability of the 
results obtained in this work. 
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